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'i^^'. Polarization observables in the *He(e, e'p)''H reaction are calculated using accurate three- and 

' four-nucleon bound-state wave functions, a realistic model for the nuclear electromagnetic current 

, operator, and a treatment of final-state-interactions with an optical potential. In contrast to earlier 

^S) ' studies, no significant discrepancies are found between theory and experiment both for the ratio of 

^ ' transverse to longitudinal polarization transfers and for the induced polarization, when free-nucleon 

^ , electromagnetic form factors are used in the current operator. The present results challenge the 

' current interpretation of the experimental data in terms of medium-modified form factors. 
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^"H ' The recent measurement, carried out at Jefferson Lab (JLab) [1] (E93-049), of the ratio of transverse, P^, to 
^ . longitudinal, P^, polarization transfer parameters in the ''He(e, e'p )'^H reaction has generated considerable interest. 

' In the elastic process ep ep, the PUP'z ratio is proportional to that of electric to magnetic form factors of the 
proton [2], and therefore its measurement in a nucleus by quasi-elastic proton knockout can shed light, in principle, on 
fS| the question of whether these form factors are modified in medium. However, the answer to this question is obviously 

, model dependent, since the modification is inferred from a comparison of the experimental data with theoretical 

■ predictions for (e, e'p) cross sections in nuclei. Therefore, it is crucial for a proper interpretation of the experimental 
, data that the theoretical calculations include contributions from final state interactions (FSI) between the knocked 

1-^5 ' out proton and residual system, as well as from many-body terms in the nuclear electromagnetic current and from 
I , correlation effects in both the initial and final bound nuclear clusters. At issue then is whether these contributions have 
Q ' so far been accounted for and reliably estimated. For example, the studies of Refs. [3,4], based on relativistic mean 
^ , field theory, ignore correlation effects in the bound state wave functions and many-body terms in the electromagnetic 
^ operator. Furthermore, FSI are treated with a relativistic optical potential in the work of Udias and collaborators [3], 
^ in which the contributions associated with charge exchange processes are neglected — they will turn out to play an 
important role in the reaction under consideration, see below. In the work of the Ghent group [4], FSI are described in 
a Glauber framework, which may not be reliable at the low end of the Q^-range covered by E93-049, since the ejected 

■ proton energies are too low. In addition and more importantly, the charge-exchange mechanism referred to earlier 
is also not included in this study — indeed, it is not obvious how to incorporate it within the context of a Glauber 
approach. Lastly, Laget's calculations [5], a full account of which is yet to be published, treat FSI by retaining S, P, and 
D waves in the nucleon-nucleon (NN) scattering amplitude at low energy, and by using a standard parameterization 
of the latter in terms of a central term at higher energies. Charge-exchange as well as spin-dependent effects beyond 
those implicit in the use of the low-energy NN amplitudes, are neglected. Two-body terms in the current operator are 
shown to lead to a quenching of ~ 2-2.5 % in the ratio P^/P^ relative to that in Plane- Wave-Impulse- Approximation 
(PWIA) — note, however, that in Ref. [1] no quenching is reported in the result at Q^=Q.5 (GeV/c)^ (the only one 
shown) for the calculation from the same author. Approximations are made in the numerical evaluations of the loop 
integrals occurring in Laget's diagrammatic approach. 

The present study is based on variational wave functions for the bound three- and four-nucleon systems, derived 
from a realistic Hamiltonian consisting of the Argonne wig two-nucleon [6] and Urbana-IX three-nuclcon [7] interactions 
(AV18/UIX) with the hyperspherical- harmonics (HH) technique, with (^He/'^H) [8] and without (^He) [9] the inclusion 
of pair correlations. The high accuracy of the HH wave functions is well documented [10], as is the quality of the 
AV18/UIX Hamiltonian in successfully and quantitatively accounting for a wide variety of three- and four-nucleon 
bound-state properties and reactions, ranging from binding energies, charge radii, and elastic form factors [10-12] to 
low-energy radiative and weak capture cross sections and polarization observables [13], to the quasi-elastic response 
in inclusive (e, e') scattering at intermediate energies [14]. 
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The polarization transfer measurement in the JLab experiment E93-049 was performed in a quasi-elastic regime: the 
momentum of the recoiling "^H nucleus was kept close to zero. The proton lab kinetic energies were (0.29,0.55,0.88,1.42) 
GeV for the values (0.5,1.0,1.6,2.6) (GeV/c)^, respectively. These energies are obviously beyond the range of 
applicability of A^A^ interaction models, such as the AV18, which are constrained to reproduce A^A'^ elastic scattering 
data up to the pion production threshold. At higher energies, A^A^ scattering becomes strongly absorptive with the 
opening of particle production channels. Indeed, the pp inelastic cross section at 0.5 GeV increases abruptly from 
about 2 mb to 30 mb, and remains essentially constant for energies up to several hundred GeV [15]. 

In view of these considerations, FSI in the p^H scattering state are described in the present work via an optical 
potential [16,17]. Of course, this approximation has limitations as to the energy range where it is expected to be 
valid, see discussion below. The p wave function is then written as 

= 5^(-)P"™[r7[;'(^;p)0.3(ifc^;'H) + r?(;H^;n)0,3(jfc^;3He) 

perm 

where <j and are the spectator nucleon and bound cluster spin projections, k is their relative momentum, and the 
sum over permutations ensures the antisymmetry of the wave function The spectator wave functions r]{i;p/n) 

are given by the linear combinations [r]{i;T = l)+/—ri{i;T = 0)]/2, where T=0,1 denotes the total isospin of the 
1+3 clusters. The latter are taken to be the scattering solutions of a Schrodinger equation containing a complex, 
energy-dependent optical potential of the form 

<*(7;ei) = [v%r; T,ei) + (4T - 3)v'=^r; T,ei)] + [v^r; T,ei) + (4T - 3)v''^r; T,ei)] I • s , 

where T^ei is the relative energy between clusters i and jkl, and I and s are the orbital and spin angular momenta 
of nucleon i, respectively. The imaginary part of accounts for the loss of flux in the p^H and n^He states due 
to their coupling to the dd, three- and four-body breakup channels of ''He. Note that the n+'^He component in 
the scattering wave function tp^~^p +^11) vanishes unless the isospin-dependent (charge-exchange) terms in v°^^ are 
included. In the results presented below, all partial waves are retained in the expansion of 77(1; T), with full accoimt 
of interaction effects in those with relative orbital angular momentum I < 17. It has been explicitly verified that the 
numerical importance of FSI in higher partial waves is negligible. 

The central v'^ and v^'^ , and spin-orbit and v^'^ terms have standard Woods-Saxon and Thomas functional forms. 
The parameters of u^, w^"^, and were determined by fitting p+^H clastic cross section data in the lab energy range 
'71ab= (160-600) MeV, and p -|-^ H ^ n -|-^ He charge-exchange cross section data at Tiab=57 McV and 156 MeV (see 
Refs. [16,17] for a listing of their values). The charge-exchange spin-orbit term is taken to be purely real, with a 
depth parameter depending logarithmically on Tjab, 15.0— 1.5 log[Tiab(MeV)] in MeV, and with radius and diffuseness 
having the values 1.2 fm and 0.15 fm, respectively. The isospin-independent and isospin-dependent spin-orbit terms 
of Uy^* are not well constrained by the data, since these consist exclusively of differential cross sections [16,18,19]. 
However, they significantly affect the induced polarization Py measured in the ^He(e, e'p )''H, as will be shown below. 

The nuclear electromagnetic current includes one- and two-body terms. The one-body current and charge operators 
have the form recently derived by Jeschonnek and Donnelly [20] (specifically. Eqs. (23) and (25) of Ref. [20]) from 
an expansion of the covariant single-nucleon current, in which only quadratic and higher order terms are neglected 
in its dependence on the initial nucleon momentum. This form of the one-body currents is well suited for dealing 
with processes in which the energy transfer may be large (i.e., the ratio of four- to three-momentum transfer {Q/q)"^ 
is not close to one) and the initial momentum of the struck nucleon is small. Thus, its use is certainly justified in 
the quasi-elastic kinematics of the E93-049 experiment under consideration. In the limit {Q/qY — 1 one recovers 
th(^ standard expressions for the impulse-approximation currents (including the spin-orbit correction to the charge 
operator) . 

The two-body charge and current operators consist of a "model-independent" part, that is constructed from the A^A'' 
interaction (the AV18 in the present case), and a "model-dependent" one, associated with the excitation of intermediate 
A resonances and the pir^ and uj-k^ transition mechanisms (for a review, see Ref. [12] and references therein). 

Finally, the Hohler parameterization [21] is used for the electromagnetic form factors of the nucleon, except at the 
highest values of 1.6 (GeV/c)^ and 2.6 (GeV/c)^, for which the proton electric and magnetic form factors are taken 
from the parameterization obtained in Ref. [22] by fitting Gup data and the ratio Gep/Gmp recently measured at 
JLab [2] . Incidentally, E93-049 also reported measurements of the polarization transfer ratio on hydrogen in the same 
kinematics as for helium [1]. The results for (P^/P^)ih are in agreement with the Hohler-Brash parameterization, 
except at (3^=1.6 (GeV/c)^: Strauch et al. measure — 0.395 ± 0.013 while the Brash fit, based on Ref. [2], gives 
—0.415. However, this difference is due to finite detector acceptances [23]. 
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The matrix elements '^^^p +"^11; k cr, (73 | j^(q, a;) |^He) are computed with Monte Carlo (MC) techniques without 
making any further approximations beyond those inherent to the treatment of FSI and nuclear electromagnetic 
currents, discussed above. The resulting theoretical predictions for the super-ratio R={P!^/P^)/{P!^/P!,)pwia and for 
the induced polarization Py arc compared with the experimental data [1] in Figs. 1 and 2. The ratio of transverse to 
longitudinal polarizations in PWIA is proportional to Gep/Gmp as obtained in the Hohler-Brash parameterization. 
In Fig. 1 the hydrogen data [23] are also shown, for which, as expected, R is very close to one, except for the point 
at Q^—l.Q (GeV/c)^ (sec comment above, however). 

The calculated results in Figs. 1-2 are labeled as follows. The curves OPT(no CH-EX) and OPT both use one- 
body currents and the optical potential to describe FSI effects, the only difference being that in the OPT(no CH-EX) 
calculation the charge-exchange components v'^'^ and v^^ of u^^* are ignored. The curve labeled OPT+MEC includes 
the full ii^f* (as the curve OPT) and one- and two-body currents. The statistical errors associated with the MC 
integrations are only shown for the OPT-I-MEC predictions, they are similar for the other predictions. Finally, the 
results of a calculation including one- and two-body currents, in which the sign of the charge-exchange spin-orbit 
term v^'^ in v':^'^ had been artificially flipped, were found to be numerically close to those obtained in the OPT 
approximation, except for Py at Q^=0.5 (GeV/c)^, i.e. Pj^=— 0.0060 ± 0.0090. They are not shown in Figs. 1 and 2. 

It should be stressed once more that the calculations may not describe reliably FSI effects for the last two 
values, since the relevant proton kinetic energies, 0.88 GeV and 1.42 GeV, represent uncontrolled extrapolations of 
the present optical model, which is fit to data up to 0.6 GeV. For the low values, however, the Py results obtained 
in the OPT approximation indicate that, while the spin-orbit terms and v^"^ may not be well constrained by the 
p^H elastic and charge exchange differential cross sections, they seem nonetheless to be quite realistic. 

The OPT-I-MEC calculation reproduces well the measured super-ratio R at low values, and is also consistent 
with the measured induced polarization Py, although the experimental data for this latter quantity have rather large 
systematic errors. The charge-exchange components v'^'^ and, particularly, v^'^ in the optical potential play a crucial 
role — see curves OPT(no CH-EX) and OPT in Fig. 1 — as do two-body terms in the electromagnetic current operator. 
The inability to reproduce the observed quenching of the super-ratio had been a persistent problem in all earlier 
studies we are aware of [3-5]. Indeed, the results of these studies are similar to those obtained here in the OPT (no 
CH-EX) calculation. 

In the parallel kinematics of E93-049 P'^ and P'^ are proportional, in the notation of Ref. [24], to the response 
functions R^i^j,, and Ripj., , respectively, involving interference between matrix elements of charge-current and current- 
current operators. In fact, it turns out that i?^y,=iZT exactly (again, in parallel kinematics), where Rt is the ordinary 
transverse response. The charge-exchange mechanism affects both R\rp, and Rlj-rp,: for example, at Q^=0.5 (GeV/c)^, 
i?^j„=0.251 fm^ and 0.227 fm^, and ii^r-^, =0.183 fm^ and 0.174 fm^ in the OPT(no CH-EX) and OPT calculations, 
respectively. However, the resulting polarization transfer parameters are: =-0.116 and -0.115, and P^=0.130 and 
0.136 in the same approximations. 

Two-body terms in the current operator (those in the charge operator give tiny contributions) also affect P'^ and 
P^ differently. Both response functions R^j-, and Rj-j.,=Rt are increased by two-body current contributions, but 
the increase for Rt, about 8 %, is twice as large as for R\rp, (as one would naively expect), and therefore the ratio 
Px/P'z is suppressed by about 4 % with respect to that obtained with one-body currents only. The enhancement of 
Rt is consistent with that calculated for the transverse response function, measured in inclusive ^He(e, e') scattering 
in quasi-elastic kinematics [14], although it is important to emphasize that the total p'^H contribution to the inclusive 
response involves an integral over the missing momentum p^, while here this contribution is evaluated at a single 
kinematical point, namely Pm — 0. Lastly, among the two-body terms, the 7r-like and p-like currents, derived from 
the isospin-dependent static part of the AV18, and the A-excitation current give the leading contributions. 

To conclude, the observed suppression of the super-ratio in ''He is explained by FSI effects and two-body current 
contributions. In contrast to earlier suggestions made in the literature [1], no in-medium modification of the proton 
electromagnetic form factors is needed to reproduce the experimental data. The present results corroborate the 
conclusions derived from analyses of the Coulomb sum rule (GSR) in few-nucleon systems [25], which show that there 
is no missing strength in the longitudinal response of these nuclei when the free-proton electric form factor is used. The 
CSR situation for medium- weight nuclei remains controversial to this day [26,27], although there are rather strong 
indications that even there no quenching of longitudinal strength is observed [26,28]. Therefore, the quark- meson 
coupling model of nucleon and nuclear structure [29], which leads to the notion of medium-modified nucleon form 
factors, seems to be at variance with a number of experimental observations. It is interesting to note that this notion 
is not an inevitable consequence of the quark substructure of the nucleon. For example, a recent study [30] of the two- 
nucleon problem in a flux-tube model of six quarks interacting via single gluon and pion exchanges suggests that the 
nucleons retain their individual identities down to very short separations, with little distortion of their substructures. 
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FIG. 1. The super-ratios measured for *He are compared to theoretical predictions, obtained in various approximation 
schemes (see text for an explanation of the notation) . The solid lines are to guide the eye only. Also shown are the super-ratios 
measured for ^H. Note that the ^He (^If) data have been shifted to the left (right) by 0.05 (GeV/c)^ in order to reduce clutter. 
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FIG. 2. The induced polarizations measured for *He are compared to theoretical predictions, obtained in various approxi- 
mation schemes (see text for an explanation of the notation) . The solid lines are to guide the eye only. Note that the ''He data 
have been shifted to the left by 0.05 (GeV/c)^ in order to reduce clutter. 
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